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Abstract —The objective of our paper is to improve efficiency 
(in terms of throughput or system capacity) for mobile satellite 
communications. In this context, we propose an enhanced Hybrid 
Automatic Repeat reQuest (HARQ) for delay tolerant services. 
Our proposal uses the estimation of the mutnal Information. We 
evalnate the performance of the proposed method for a land 
mobile satellite channel by means of simulations. Results are 
compared with those obtained with a classical Incremental re¬ 
dundancy (IR) HARQ scheme. The technique we propose, shows 
a better performance in terms of efficiency while maintaining an 
acceptable delay for services. 

Keywords—Hybrid ARQ; Satellite Communications; Land Mo¬ 
bile Satellite (IMS) Channel; Delay; Efficiency. 

I. Introduction 

For many years, mobile satellite communications services 
are challenging and very costly. Land mobile satellite (LMS) 
channels are highly affected by important propagation impair¬ 
ments (time-selective channel) and inter/intra system interfer¬ 
ence, that both cause unstable and low signal to noise ratio 
(SNR), which the receivers have to cope with. This usually 
implies strong limitations on the delivered service throughput, 
whatever the considered satellite frequency band. 

Our objective is to propose a mechanism, which improves 
the efficiency of link usage while providing an appropriate 
service to applications. The targeted services (data transfer 
from sensors, messages for aeronautical services, etc.) are 
assumed to be tolerant to delay. For example some aeronautical 
services define delay requirement for the delivery of 95% of 
messages |[T). 

To deal with problems caused by link characteristics in 
mobile satellite communications, there are many solutions. 
One of these solutions is to use pure Forward Error Correction 
(FEC), which can makes the message very robust. However 
using FEC is not sufficient sometimes, due to the highly 
varying channel, where its quality changes dramatically. This 
make it difficult sometimes to decode the message, even if 
the used code is very robust. Alternatively, Automatic Repeat 
reQuest (ARQ) can be used as a solution to deal with channel 
variations, where transmitter attempts many retransmissions in 
case of unsuccessful decoding. ARQ messages without FEC, 
are not so robust to ensure a reliable communication. 


Hybrid ARQ (HARQ) protocols are used in most of re¬ 
cent terrestrial wireless communication systems. Worldwide 
interoperability for microwave access (WIMAX) and long 
term evolution (LTE) are examples of these systems that 
use HARQ Q. There are many types of HARQ techniques. 
HARQ type I is called chase combining (CC). The receiver 
asks the transmitter for the transmission of the same packet 
of coded data. At the receiver, the decoder combines the 
multiple copies of the transmitted packet weighted by the 
SNR received. HARQ type II, is also called incremental 
redundancy (IR). Unlike the previous method, which transmits 
simple repetitions of the same encoded packets, IR technique 
transmits additional redundant information in an incremental 
way if the decoding does not succeed from the first trans¬ 
mission. If each retransmission packet is self-decodable this 
scheme is called type III HARQ ||^. Classical IR HARQ, 
which transmits a fixed number of bits at each transmission, 
is not optimal from the efficiency point of view even if it can 
improve the decoding delay especially for high values of SNR. 
However, choosing optimal values for number of bits to be 
sent at each transmission can improve the system performance 
and throughput level. Many papers has proposed methods to 
choose an optimal number of bits at each transmission, a 
theoretical basis for a parity bits selection by means of a risk- 
sensitive optimal control is established in 0 - The idea of 
random transmission assignments of the mother code bits was 
introduced in 0 - An optimized IR HARQ schemes based on 
punctured LDPC codes over the EEC was proposed in 0 - 

Contrarily to classical IR HARQ scheme, the main idea 
of our enhanced HARQ technique is to estimate the average 
number of bits to be transmitted at each transmission to decode 
the codeword with a targeted probability. This probability 
depends on the application/service (delay constraints). This 
technique uses the mutual information to predict the mean 
number of bits needed for each transmission. It uses the knowl¬ 
edge of the statistical distribution of the channel attenuation. 
Our enhanced HARQ transmission proposal is simulated in a 
satellite communications environment, where an LMS channel 
and a long Round-Trip time are considered. 

The remainder of this paper is organised as follows. Section 
II describes the channel capacity and the LMS channel model. 


Data bits 


We present classical IR HARQ in Section III. In Section IV, 
we propose an enhanced HARQ model for delay tolerant 
services in mobile communications. We present the results 
of simulations and we compare the performance of both 
techniques (classical IR and enhanced HARQ) in Section V. 
We conclude our study in Section VI. 

H. Channel modelling 
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0bit to be sent at the Ist transmission 9^')^ be sent at the 2nd transmission 

^bit to be sent at the 3rd transmission sent at the 4th transmission 


Figure 1. Example of the transmission technique considered in our 
simulations (Code rate 1/6, with a maximum of 4 retransmissions). 
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In our study on mobile satellite communications, we con¬ 
sidered an LMS channel to model this environment 0 ©■ In 
the following, we explain how to compute the capacity of this 
channel. 

A. Channel capacity and Mutual Information 

Channel capacity C quantifies the maximum achievable 
transmission rate of a system communicating over a band- 
limited channel, while maintaining an arbitrarily low error 
probability. It corresponds to the maximum of the mutual in¬ 
formation between the input and output of the channel, where 
the maximization is done with respect to the input distribution. 
Mutual information (MI) measures the information that input 
of the channel (X) and output of the channel (V) share. It is a 
key parameter in our approach, as it will be used to calculate 
the mean number of bits needed to decode a message at each 
transmission with a given probability. 

Given the channel input symbol Xi, its energy Esi, a 
realisation of noise rii (which has a Gaussian distribution 
with variance and the channel attenuation pi, the channel 
output symbol pi can be written as; 

yi = Pi\/~E7iXt + rii. ( 1 ) 

For an equally distributed input probability, the MI corre¬ 
sponds to the capacity of the channel, which is AWGN, can 
be calculated by the following equation Q: 
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where; 

• M is the modulation order; 


• is the space dimension that depends on the used 
modulation (N = 2 for any PSK-based modulation with 
more than 2 states); 

• f is the integration variable of dimension N (t = 


For the rest of the paper, we define Ml^eq as the average MI 
per bit required to decode a codeword at a given probability 
expressed in Word Error Rate (WER). MIreq can be calculated 
using the mutual information function Q and the performance 
curves giving the WER versus Es/Nq of the used modula¬ 
tion/coding scheme m- By using these performance curves, 
we can deduce the Es/Nq necessary to obtain a given WER. 
Then we use this value of Es/Nq in the mutual information 
function 0 to obtain the mutual information required to 
decode a codeword at this given WER. The prediction of 
performance of the WER based on MI is quite classical, and 
has been described and validated in GD- 

Erom now on, to compute the mutual information we use 
0, ^ is replaced by Pi^, where ^ is the average energy 
per symbol. 

B. LMS Channel 

One of the reference propagation models for LMS channel 
is a statistical model based on a three state Markov chain Q. 
This model considers that the received signal originates from 
the sum of two components; the direct signal and the diffuse 
multipath. The direct signal is assumed to be log-normally 
distributed with mean a (decibel relative to LOS (Line Of 
Sight)) and standard deviation T' (dB), while the multipath 
component follows a Rayleigh distribution characterized by its 
average power, MP (decibel relative to LOS). This model is 
called Loo distribution Eor the modelling of the LMS 

channel in our simulations, we use attenuation time series 
using a propagation simulator based on the three state channel 
0 0 provided by CNES. Using this tool we calculate the 
distribution of the probability to obtain an attenuation in the 
channel for a given environment. 

HI. Incremental Redundancy HARQ 

IR HARQ can be described as follows. The sender transmits 
a number of bits that correspond to a given codeword. After 
receiving the feedback (ACK/NACK) from the receiver, the 
transmitter decides to no longer send bits corresponding to 
this codeword if an ACK is received, or to send more parity 
bits if a NACK is received. The number of bits to be sent in 
the next retransmission is fixed at each transmission. These 
numbers of bits are chosen without any knowledge about 
the channel quality or the global channel characteristics. If 
the cumulative received sequence can not be decoded after 
a maximum number of bits transmitted, the transmitter stops 
sending bits that correspond to this codeword. The parity bits 













are generated according to a coding scheme with a code rate 
corresponding to the maximum number of bits that can be 
transmitted per codeword. Bits to be sent at each transmission 
are part of the original codeword (mother code), leading to a 
different code rate at each transmission (see Figure [T]). This 
technique of transmission is somehow similar to puncturing. 
Data bits of each codeword must be kept at the buffer of 
the sender as long as the codeword is still not decoded or 
the transmitter has not decided to end the transmission. At 
the receiver side, the LLRs (Log Likelihood Ratio) of the 
received symbols of each codeword are kept in the buffer as 
long as the codeword is not decoded or the transmission of 
the corresponding bits has not ended. An overview of ARQ 
and HARQ mechanisms implemented or proposed in beyond 
3rd generation systems was presented in Q. 

IV. Enhanced HARQ for delay tolerant services 

IN MOBILE SATELLITE COMMUNICATIONS 

In this section, we propose an enhanced HARQ model using 
mutual information. In our scheme, the computation of the 
numbers of bits to be transmitted at each transmission is done 
in a way that insures the decoding with a probability. These 
probabilities are predefined for each transmission and chosen 
according to delay constraints of the application. 

In the following, we explain how to proceed to compute 
the number of bits required to decode a codeword with a 
predefined probability, supposing the knowledge of the global 
statistics of the channel. 

A. Enhanced HARQ Model 

Our proposal uses the MI to compute the number of bits to 
be transmitted at each transmission. After transmitting some 
bits of a given codeword, this codeword accumulates some 
mutual information. This mutual information can be computed 
knowing the number of bits transmitted and the attenuation 
coefficient affecting the transmitted bits. In this model we 
consider a reference Es/Nq, which is a fixed value of Es/Nq 
in clear sky (for a given terminal and without attenuation). 
The model assumes that the channel is stationary for the 
transmission time of the bits, at a given transmission for a 
given codeword. The MI obtained at the transmission for 
a given codeword can be computed as: 




where: 

• is the attenuation coefficient affecting bits transmitted 
at the transmission for a given codeword; 

• ^serit '^he number of bits transmitted at the trans¬ 
mission; 

• MI{-) is the function giving the value of mutual infor¬ 
mation for a given Es/Nq on a gaussian channel. 


The MI per bit accumulated for a given codeword, from the 
beginning of transmission until the transmission, can be 
computed as: 
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where: 

. Af(j) is the total number of bits transmitted for a code¬ 
word up to the transmission; 

Note that is unknown in our model, thus MlaJc is also 
unknown. 

Let us consider the minimum MI per bit needed 

to decode the codeword at the (j -f 1)*^ transmission with a 
predefined decoding probability. 

Let be the number of bits to be transmitted at the 

(j + 1)*^ transmission, we have: 

N,usMNs, = , ( 5 ) 


where N^its is the maximum number of bits that can be 
transmitted for a codeword in total. 

Finally is given by: 

,^(, + 1 ) _ NutsMIreq - 

needed S/rrO + l) ' 

^needed 

^^nexded ^laX are the key parameters for the com¬ 

putation of the number of bits to be sent at the (j + 1)*^ 
transmission. In the following, we explain in detail how to 
proceed to calculate these two values at each transmission 
according to the predefined decoding probabilities. 

B. Computation of and MlaX 

In a first step, we will explain how to calculate 
then we finish by explaining the way we compute 

We use the knowledge of the statistical distribution of the 
channel attenuation to control the probability of decoding a 
codeword at each transmission. 

The idea is to define at the beginning of the communica¬ 
tion, a table containing the probability of decoding at each 
transmission. The decoding probability and the efficiency of 
the link usage are related, the sender can transmit a large 
number of bits at the first transmission, which increases the 
decoding probability but the efficiency will decrease and vice 
versa. So, we have to improve efficiency while respecting 
delay constraints for services. In addition the sender may also 
want to limit the number of transmission attempts for a given 
codeword. 

In the rest of this paper, we will consider Pj the probability 
of decoding at the transmission conditioned on the fact 
that decoding at earlier transmissions was impossible, where 
P = Pj is the percentage of decoded codewords over all 
the transmitted codewords and is equal to 1. 
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To target a decoding probability Pj at the transmission, 
we have to find the corresponding necessary to 

calculate the number of bits to be transmitted 

For this, we define as the minimum successive atten¬ 

uation coefficient (threshold) that provides a Ml greater than 
Ml^^eded- threshold is the minimum attenuation needed 
to obtain the decoding probability at the transmission. 
Pneeded d^pends not Only on the element in the probability 
decoding table but also on Pk as we will see later in 

^Inleded given by: 

Mllilded = ( 7 ) 

Where MI{-) is defined in 'll, that takes as input 

( U) ,12 

\r needed' 'No' 

We assume that the reference Es/Nq fixed for a given ter¬ 
minal (only change over the time). depends only 

on the channel attenuation threshold Since decoding 

probabilities at each transmission are predefined, p^^leded i^ 
pre-computed. MI{-) is a strictly increasing function (as a 
function of p). Then any attenuation coefficient greater than 
Pnleded vviH lead to a successful decoding. To determine 
Pneeded leading to Pj, we use the cumulative distribution 
function (CDF) of the attenuations of LMS Channel. 

To simplify our calculation, we consider these two events: 

• Aj: Successful decoding at the transmission; 

• Bj-i'. Not decoding at the {j — 1)*^ transmission. 

Pj can be defined as p{Aj n Sj-i) and pj is p{Aj). Since 
Aj and Bj-i are independent (according to the channel 
modelling), 

fc=i (8) 

= 

CDF of the channel gives us P{p < p^^leded)^ while 
Pj corresponds to P{p > (successful decoding). 
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= 11010 bits 


Using (6) 


Mrall = 'ZiP„i MiapiY = 1.4104 / (Pi < -4.1861) Using (9) 
^0 




^bits 21761-11010*1.4104 


bits 

MlS^c = I.iPwiMI((Piy ^) = 0.8353 / (Pi <-5.9094) 
^0 




Nbits 21761-14278*0.8353 


= 33623 bits 


Figure 3. Numerical example for the computation of number of bits to be 
transmitted at each transmission (Es/Nq = 13 dB, code (8920,1/6), three 
transmissions). 
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Therefore, p^^eded CDF graph is given by l-pj (See 

Figure |^. 

Once we found pnleded l^^ding to Pj, we use it in <0 to 
calculate 

At the (j — 1)*^ transmission, only codewords affected 
by attenuation coefficients greater than are decoded. 

AIIaic^'^ represents the average mutual information per bit of 
all codewords that have not been decoded at the (j — 1)*^ 
transmission, affected by attenuation coefficients less than 
Pneeded- ^^acc is the sum of all mutual informations 
weighted by their probabilities in the non decoding zone (see 
Figure]^. For this, we use the CDF and the PDF (probability 
density function) of the channel. Mlacc can be calculated 
by: 

MNiiN = ( 9 ) 

/Vq 
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where: 

• Pi belongs to the ensemble of attenuation coefficients, 
which are less than p^^di 

• Pwi is the weighted probability, defined as the probability 
to obtain pi in the non decoding zone. 

Finally, to calculate we use and MliiN 

in (|^. Where iV'-^ is the sum of all numbers of bits already 
calculated for the previous transmissions. 

Figure represents a numerical example for computation 
of the numbers of bits to be transmitted at each transmission, 
showing different steps for the calculation of all parameters in 
our model. 


Figure 2. CDF of the attenuation coefficients in the LMS channel. 












































V. Simulations and Results 

The system parameters that we will consider for all simu¬ 
lations are; 

• Satellite Orbit; Geostationary (GEO) 

• Round Trip Time; 500 ms 

• Band; S 

• Land mobile satellite channel , Intermediate Tree Shad¬ 
owed Environment (ITS) 

• Speed; 60 Km/h 

• Distance; 10 Km 

• Mother EEC code, CCSDS Turbo Codes 1/6 

• Codeword length; 53520 bits (Data bits ; 8920 bits) 

• Modulation; QPSK 

• Symbol time ; 4.10“® seconds, bit rate (Rb); 500 Kbps 

Our simulations are about 10 minutes of communication 
between the transmitter and the receiver (about 300 Mb 
transmitted). In our simulations, we consider a targeted WER 
of 10“^, and we use the actual performances of CCSDS 
Turbo codes (8920,|) as presented in p^ . Attenuations in ITS 
environment are very high, which makes difficult the decoding 
of 100% of codewords in the range of reference Es/Nq 
considered (7 to 13 dB) with code rates greater than |. Hence 
the choice of code rate g. We suppose that the synchronisation 
is never lost. The transmitted codewords are identified with 
a sequence number that is never lost. We suppose also that 
the return channel does not introduce errors and the feedback 
can be transmitted immediately (no congestion problem on the 
reverse link). 

A calibration phase is required in the simulations, that take 
into account the actual numerical performances of the targeted 
EEC code(s). This allows to avoid implementing a real decoder 
in the simulation chain, while assuring a very good accuracy 
of the representation Ol- We use the MI to decide if the 
codeword is decoded or not using this formula; 

iV0)M/« > NkusMIreg . (10) 

Simulations of classical IR and proposed enhanced HARQ 

We present results obtained by implementing both schemes 
(classical IR and proposed enhanced HARQ) described in 
previous sections, and we compare these results. 

We define a set of precomputed decoding probabilities, that 
provide for each retransmission a precomputed number of 
parity bits to send. We are interested to improve the efficiency 
while maintaining an acceptable delay for services. 

We define the efficiency, E (bits/symbol), as follows; 

^data bits^decoded words /ii\ 


where; 

• Ndata bits Is the number of data bits (useful bits) per 
codeword considered in our coding scheme; 


• -^decoded words Is the total number of decoded codewords 
during the communication; 

• Ntotai is the total number of symbols transmitted during 
the communication. 

The delay for decoded codewords (at the receiver) can be 
expressed in terms of number of transmissions (Ntrans), bit 
rate (Rb), number of bits sent (N) and propagation delay 
(Tpropag), assuming a negligible access delay; 

N 

Delay = —— -t- ‘2,(JS[tj-ans l)^^ropa5 Tpropag (-S) . (12) 
Kb 

The efficiency and the delay are positively related. To 
improve the efficiency we can increase the delay while 
respecting the delay constraints for services. This delay is 
controlled by the predefined decoding probabilities, which will 
be quasi constant along all the values of reference Es/Nq. In 
the following, we have considered three different tables of 
decoding probabilities. Eixed decoding probabilities at each 
transmission considered in our simulations are given in Table 

where the maximum number of transmissions for a given 
codeword is four. Case 1 corresponds to a service accepting the 
delivery of 80% of the messages at the first two transmissions 
and 20% at the last two retransmissions. Case 2 corresponds to 
a service accepting the delivery of 95% of the messages at the 
first two transmissions and 5% at the last two retransmissions. 
Case 3 corresponds to a service accepting the delivery of 45% 
of the messages at the first two transmissions and 55% at the 
last two retransmissions. 

TABLE 1. PREDEEINED DECODING PROBABILITY TABLE FOR THE 
3 CONSIDERED CASES 


Transmission 



ya 

^tfl 

Pi{casel) 

0.5 

0.3 

0.15 

0.04999 

Pi(case2) 

0.75 

0.2 

0.03 

0.01999 

Pi{case3) 

0.4 

0.05 

0.45 

0.0999 


The number of bits to be transmitted at each transmission for 
the classical IR HARQ scheme (Section III) can follow several 
strategies. As a simple case, we consider an equally shared 
repartition of the data-nparity bits among the transmissions, as 
shown in Table [II] 

TABLE II. NUMBER OE BITS TO BE TRANSMITTED AT EACH 
TRANSMISSION FOR CLASSICAL IR HARQ MODEL 


Transmission 


<2nd 

yd 

^th 


13380 

13380 

13380 

13380 


After simulating first (classical IR HARQ) and second (en¬ 
hanced approach with three cases mentioned above) schemes, 
under the same conditions, we figure out the efficiency and 
the delay obtained for both schemes and compare them. 

Eigure 1^ shows the average delay, required to decode code¬ 
words, obtained with both models (by means of simulations). 
The mean delay is computed by averaging delays obtained 
for decoded codewords calculated using As we can see, 
the values obtained for each case of the proposed model 
are approximately the same. The delay remains stable and 



















Figure 4. Delay obtained for both models for different values of reference 

iEs/No)dB- 

controlled by the decoding probability, globally constant. It 
changes a little, from a reference Es/Nq to another, according 
to the number of bits transmitted and the values of decoding 
probabilities. 

Figure [^compares the efficiency obtained with both models 
in the same conditions. We can see that even if we loose 
in delay, the proposed model in cases 1 and 2, outperforms 
the classical one especially for high values of reference 
Es/Nq. The resulting gain can reach 15%. These results 
validate our proposal, that fixing decoding probabilities at each 
transmission can improve the efficiency while respecting delay 
constraints for a given service. These results seem promising, 
however they have been obtained without any optimization 
of decoding probability for each transmission step. Therefore, 
some further improvements could probably be obtained. This 
clearly calls for an optimization process in further steps of the 
work. 

VI. Conclusion AND Future Work 

In this paper, we have compared two techniques of HARQ 
transmission. The hrst one is a classical IR HARQ scheme, 
that transmits the same number of bits at each transmission 
without any knowledge about the global statistics of the chan¬ 
nel; the second one is an enhanced HARQ technique, which 
takes into account the global statistics of the channel. This 
proposed enhanced technique estimates the number of bits to 
be transmitted at each transmission to decode a codeword with 
a given probability. This estimation relies on the mutual infor¬ 
mation given the decoding probability at each transmission and 
the knowledge of the distribution of the channel attenuations. 
Finally, results obtained after simulating both schemes in a 
mobile satellite communication environment are compared in 
terms of decoding probability (delay) and efficiency. Results 
show that enhanced HARQ has better performance in terms 
of efficiency especially for high values of reference SNR, 
while maintaining a quasi constant delay acceptable by delay 
constraints of service. As future work we will modify the 
proposed method to make it an adaptive model. This model 



E„/N„(dB) 

Figure 5. Efficiency obtained for both models for different values of 
reference {Es/No)dB- 

would calculate the number of bits to be transmitted at 
each transmission according to decoding probabilities and 
real values of attenuation coefficients obtained. We plan also 
to consider additional system parameters, such as framing 
constraints and overhead to evaluate the performance of the 
mechanism in real systems. 
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